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1 .  Introduction 

The  determination  of  the  line  of  separation  over  three-dimensional 
bodies  of  arbitrary  shapes  is  a  difficult  problem  requiring  the  solution 
of  the  boundary-layer  equation  or  the  full  Navier-Stokes  equations.  A 
good  number  of  careful  boundary-layer  calculations  have  been  carried  out 
in  the  past  decade  for  the  flow  about  a  prolate  spheroid.*-5  These  cal¬ 
culations  have  clearly  demonstrated  that  the  difficulty  of  the  problem 
centers  mostly  about  the  proper  choice  of  the  grid  shape.  Most 
recently,  further  improvements  on  methods  of  integration  that  conform 
with  the  principles  of  influence  and  dependence  have  appeared. The 
difficulties  multiply  as  one  approaches  the  line  of  separation.  Taking 
into  account  the  errors  involved  in  solving  numerically  the  boundary- 
layer  equations  in  the  vicinity  of  separation  and  the  fact  that  these 
equations  are  not  even  valid  in  this  region,  the  usefulness  of  such 
methods  becomes  questionable.  Moreover,  their  application  to  more 
general  body  shapes  is  a  formidable  exercise,  requiring  the  constant 
attention  of  experienced  investigators. 

One  of  the  purposes  of  the  present  effort  is  to  provide  or  expand 
upon  some  significantly  more  simple  methods  of  calculation  of  the  separ¬ 
ation  line.  The  first  is  based  on  local  similarity  concepts  and  permits 
the  estimate  of  the  skin  friction  distribution  as  well  as  the  line  of 
separation.  The  second  allows  the  determination  of  the  separation  line 
by  criteria  based  exclusively  on  potential  flow.  This  is  admittedly  a 
very  simplistic  point  of  view,  but  the  accuracy  of  the  results  has  sur¬ 
prised  the  authors  and  certainly  deserves  to  be  communicated  with  the 
engineering  community. 


Q 

Approximate  methods  are  available,  as  for  example  integral  methods 

q  in 

or  streamline  methods.  ♦  These  methods  are  considerably  simpler  than 
the  methods  based  on  the  complete  three-dimensional  boundary- layer  equa¬ 
tions.  The  new  methods  proposed  here  are  even  simpler  and  could  easily 
be  programmed  on  a  pocket  computer  to  provide  a  quick  and  reasonably 
accurate  answer  for  a  complex  configuration.  However,  a  more  promising 
approach  would  be  to  couple  a  streamline  method  *  with  a  good  local 
similarity  model.  In  this  way  the  history  of  the  developing  boundary 
layer  can  be  accounted  for.  In  this  report  a  first  attempt  along  this 
direction  is  included. 

The  second  part  of  this  report  includes  some  flow  visualizations  of 
the  same  problem.  Although  earlier  experiments  have  been  reported, 
this  experiment  was  conducted  for  Reynolds  numbers  below  2000,  which 
ensures  that  the  flow  is  purely  laminar  and  free  of  any  growing  insta¬ 
bilities.  Comparison  of  experimental  data  with  the  theory  are  included. 
However,  a  careful  comparison  will  require  digitizing  the  flow  visuali¬ 
zation  information.  The  experimental  results  offer  some  very  exciting 
new  information  on  the  physical  properties  of  three-dimensional  separa¬ 
tion  and  allow  the  numerical  results  obtained  by  others  to  be  re¬ 
interpreted. 

In  the  last  section  of  the  report,  a  two-  and  three-component  LDV 
rig  especially  designed  and  constructed  for  the  measurement  of  three- 
dimensional  boundary  layers  is  described.  Experimental  data  are  pre¬ 
sented  for  two-component  velocity  profiles  at  several  selected  stations. 
These  data  are  compared  to  boundary  layer  calculations  reported  in  Ref. 
20.  It  should  be  emphasized  here  that  to  the  knowledge  of  the  authors, 
no  other  velocity  profile  measurements  over  a  prolate  spheroid  at  such 


low  Reynolds  numbers  have  been  reported  so  far.  In  fact,  the  only 
experimental  information  on  velocity  distributions  on  a  prolate  spheroid 
pertains  to  the  cross  section  x/L  =  0.64  only,  as  reported  in  Ref.  14. 
Detailed  information  on  three-dimensional  velocity  profiles  were 

Ol 

obtained  by  Ramaprian  et  al.  by  3-hole  probes,  but  on  an  arbitrarily 
defined  body  of  revolution.  Moreover,  both  studies14,21  were  conducted 
in  the  turbulent  regime.  The  Laser-Doppler  Velocimeter  is  employed  here 
for  the  first  time  to  measure  three-dimensional  boundary-layer  profiles. 
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2.  Three  Simple  Methods  of  Calculation 

In  this  section  we  describe  three  very  simple  methods  for  estima¬ 
ting  the  approximate  form  of  skin  friction  lines  and  lines  of 
separation.  The  first  is  based  on  an  assumption  of  local  three- 
dimensional  boundary-layer  similarity.  The  second  is  based  on  some 
observations  of  existing  data.  Neither  of  the  two  requires  integration 
of  differential  equations,  although  the  first  could  be  expanded  to 
include  the  history  of  the  boundary  layer  and  integration  along  a 
streamline.  The  third  is  a  refinement  of  an  existing  streamline  method. 

For  a  few  very  special  three-dimensional  configurations  it  is  rela¬ 
tively  easy  to  obtain  accurate  numerical  solutions  of  the  three-dimen¬ 
sional  boundary- layer  equations.  One  class  of  such  configurations  are 
cylinders  of  arbitrary  cross  section  in  yaw  with  respect  to  the  oncoming 
flow.  For  these  special  bodies,  referred  to  here  as  Sears  bodies,  the 
two  momentum  equations  are  uncoupled  and  the  solution  can  be  easily 
obtained.16  We  assume  that  the  flow  about  a  general  3-D  configuration 
is  locally  similar  to  the  flow  about  one  of  the  special  bodies  discussed 
above.  In  other  words,  we  assume  that  the  properties  of  the  boundary 
layer  at  a  specific  station  on  the  body  of  interest  are  similar  to  the 
properties  at  a  suitably  chosen  point  on  a  Sears  body.  The  criterion 
for  the  determination  of  the  point  on  the  special  body  which  corresponds 
to  the  point  of  our  interest  is  based  on  properties  of  the  outer  flow. 
The  known  special  solution  is  thus  considered  essentially  a  bank  of 
properties  which  are  drawn  according  to  certain  rules,  in  order  to 
synthesize  the  properties  of  the  flow  about  the  body  of  interest.  One 
of  the  physical  and  mathematical  properties  of  the  boundary  layer  which 
is  missing  from  such  a  model  is  history.  However,  the  limitations  are 


.  S, 


i'VaV.V 


5 


not  so  severe  if  there  are  no  sharp  changes  in  the  pressure  gradient 
distribution. 

Consider  the  flow  about  a  cylinder  in  yaw.  Let  the  x,  y  and  z  axes 
be  parallel  to  the  generators,  along  the  contour  of  the  body  but  perpen¬ 
dicular  to  the  generators,  and  normal  to  the  body,  respectively  (see 
Fig.  la).  Let  u,  v  and  w  be  the  corresponding  velocity  components, 
while  U  and  V  are  the  corresponding  outer  flow  velocity  components  and 
U  ,  V<o  are  the  components  of  the  undisturbed  free  stream,  respectively. 
The  governing  equations  in  their  dimensionless  stretched  form  then  read 


a  v  +  5w  _  Q 
dy  5z 


du  .  5u  52u 

v  +  w 

ay  az  *,2 


v5-¥.+  w  —  =  V  —  + 

ay  az  dy 


In  these  equations  the  quantities  w  and  z  contain  the  stretching  factor, 
namely,  the  square  root  of  the  Reynolds  number.  Notice  that  there  is  no 
pressure  gradient  term  in  Eq.  (2).  Moreover,  Eqs.  (1)  and  (3)  can  be 
solved  independently  of  Eq.  (2).  The  familiar  boundary  conditions  are 


u  =  v  =  w  =  0 


at  z  =  0 


u  U  ,  v  -*•  V 


as  z  -*•  ® 


Solutions  to  this  problem  were  obtained  by  Sears10  for  an  outer  flow 


given  by 


v  =  Ciy  +  C3y3 


(7) 


The  solution  was  obtained  in  the  form  of  expansions 


u  «  Ujgofz)  +  jX  y2g2(z)] 


v  *  (^y  f i ( z )  +  C^3  f3(z) 


(8) 

(9) 


and  the  functions  f i ,  f3 ,  go  and  g2  were  tabulated. 

The  rules  of  correspondence  between  points  on  the  body  of  interest 
and  a  Sears  body  are  the  following: 

a.  The  distances  s  along  inviscid  streamlines  scaled  with  the  to¬ 
tal  length  of  each  streamline  are  matched, 

b.  the  U  velocity  component  is  matched, 

c.  the  V  velocity  component  is  matched,  and 

d.  the  velocity  gradient  along  the  inviscid  streamline  is  matched 
For  a  coordinate  system  x,$  defined  on  a  prolate  spheroid  as  shown 

in  Fig.  lb,  the  four  conditions  read 

x  =  x$  (10) 

v(x,$)  =  Cxxs  +  C3x$3  (11) 

u(x,*)  =  V  (12) 

<iJ-  ■  C.  ♦  3C,»sa  (13) 


where  x  is  the  reduced  distance  along  the  Sears  body. 


Equations  ( 10) - ( 13 )  are  used  to  define  a  Sears  body  and  determine  a 


poiru  on  this  body.  The  solution  at  this  point  is  then  given  by  Eqs. 
(8)  and  (9)  and  is  assumed  to  hold  at  the  corresponding  point  on  the 
body  of  interest. 

The  potential  flow  about  a  prolate  spheroid  at  an  angle  of  attack 
is  known  in  a  closed  form  solution^. 


u  _ _ 1 _ 

V®  ^_e2^X  _  t  ^2  jl/2  {(]  +  ka^cos  a>  * 

[2I_  (l)2]l/2  +  (|)(l.Kc)(Sin  a)(f  -  l)cos<>} 


(14) 


y-  ■  (1  +  Kc)sin  a  sin  *  (15) 


where  e  is  the  eccentricity 

e  =  /l-(b/a)2  (16) 

and  ka  and  kc  are  the  following  parameters 

ka  (ye  10  T^e  '  1 '  7e  **  T^f)  (17) 

i  -e 

kc  *  T+5K7  (18^ 

fl 

The  same  method  was  employed  for  the  calculation  of  the  viscous 
flow  over  an  ellipsoid  with  arbitrary  axes  a,b,c.  In  terms  of  coordin¬ 
ates  X,Y  and  Z  aligned  with  the  three  major  axes  a,b  and  c  of  the  ellip¬ 
soid,  the  potential  velocity  on  the  surface  of  the  body  can  be  calcu¬ 
lated1^  in  terms  of  an  ellipsoid  coordinate  \  which  is  constant  on  all 


confocal  ellipsoids  and  zero  on  the  ellipsoid  of  Interest.  The  solution 
is  given  in  terms  of  three  Integrals 


K-r7rr:  .fiWA'  S'A  V-  ’W 


8 


00 


d X 

— 2 - 

(a  +x  )A 


dx 

— 2 - 

(b  +X)A 


I3  -/ 


dx 

~T 


0  (a  +x)a 


(19) 


with 


A2  =  (a2+X)(b2+x)(c2+X)  (20) 

The  constants  ao.Po  and  yo  and  A,B  and  C  are  also  defined  as  follows 


oo  =  abcli , 

A  -J- 
2-ao 

(21) 

Po  a  abcl2. 

B  =  2 

2“Po 

(22) 

yo  *  abcl3. 

C  =  2 

2-yo 

(23) 

The  velocity  on  the  body  surface  is  then  given  in  terms  of  the  compon¬ 
ents  of  the  oncoming  stream,  (U  ,V  ,W  )  and  the  unit  vector  normal  to 
the  ellipsoid  at  the  point  of  interest  n^  with  components  (*,m,n) 

q  «  (AU  ,BV  ,CW  )  -  (AU  A  +  BV  m  +  CW  n)n  (24) 

The  unit  vector,  n",  normal  to  the  ellipsoid  at  the  point  X,Y,Z  has 
the  cartesian  components 


x  -  xa-2(4+4+4rl/2 

(25) 

a  b  c 

m.yb-*(4+4*4)-,/2 

(26) 

a  b  c 

7  -2,X2  Y2  Z2>-l/2 
"  ’  Zc  ("I  +  -T  *  -TT> 

(27) 

a  b  c 

In  terms  of  the  quantities,  the  cartesian  components 

U,V,W  of  the  velo- 

city  q  (see  Fig.  1  for  notation)  are  given  by 

U  •  AU  -(AU  *2  +  BV  mt  ♦  CW  rut) 
00  00  00  00 


(28) 
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(29) 

(30) 


•dd 


V  =  BV  - (AU  im  +  BV  nT  +  CW  run) 

W  =  CW  -(AU  in  +  BV  mn  +  CW  n2) 

OO  00  OO  00 

To  obtain  the  components  u,v  and  w  along  the  system  x,$,z»  one  needs,  of 
course,  to  go  through  the  proper  transformation. 

A  computer  program  was  prepared  for  the  calculation  of  viscous  and 
inviscid  streamlines.  Marching  was  started  at  points  equally  spaced  on 
a  meridional  line  and  was  marched  in  both  the  upstream  and  the  down¬ 
stream  direction.  In  this  way  difficulties  encountered  in  the  neighbor¬ 
hood  of  stagnation  points  were  avoided.  In  fact,  the  method  provides  a 
clear  indication  of  accuracy  at  stagnation  points  where  it  was  demon¬ 
strated  that  all  lines  pass  within  a  distance  of  a/1000  of  each  other. 

Sample  results  are  shown  in  Figs.  2,  3  and  4  for  a  prolate  spheroid 
and  an  ellipsoid,  respectively.  In  these  figures  are  plotted  the 
projection  of  lines  on  planes  that  correspond  to  leeward  views  and  side 
views.  The  lines  that  are  behind  the  body  are  shown  dotted.  In  all 
figures,  envelopes  of  skin  friction  lines  are  clearly  developed,  thus 
defining  the  line  of  separation.  The  separation  lines  found  are 

3 

compatible  with  analytical  predictions  of  open  separation,  however,  the 
present  method  is  very  crude  for  offering  verification  of  the  existence 
of  open  or  closed  separation.  Comparison  of  the  relative  performance  of 
this  method  with  other  methods  will  be  presented  later  in  this  report, 
together  with  the  second  method  of  separation  estimates. 

A  more  correct  approach,  which  includes  the  history  of  a  developing 
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boundary  layer,  is  based  on  equations  along  the  viscous  streamlines  * 


TO 
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(32) 

=  cosa  (33) 

Here  0  is  the  slope  of  the  streamline  with  respect  to  the  meridional 

line,  s  is  the  distance  along  the  streamline,  x„  and  x  .  are  the  shear 

x 

stress  components  along  the  x  and  the  $  direction,  respectively;  q  is 
the  magnitude  of  the  velocity,  and  f  is  the  radius  of  the  body  at  the 
location  under  consideration. 

This  system  of  equations  can  be  integrated  along  a  viscous  stream¬ 
line,  at  a  certain  elevation  from  the  wall,  provided  that  an  appropriate 
model  for  the  shear  stresses  is  available.9*10  The  shape  of  viscous 
streamlines  can  thus  be  determined.  It  would  be  desirable  to  integrate 
on  the  wall  in  order  to  obtain  the  skin  friction  lines.  Unfortunately, 
the  right-hand  side  of  Eq.  (31)  becomes  indeterminate  at  z=0,  since  both 
the  numerator  and  the  denominator  tend  to  zero.  To  circumvent  this  dif¬ 
ficulty,  a  comparison  was  made  with  a  locally  similar  flow  along  a  Sears 
body  at  each  point  of  the  numerical  integration.  In  effect  then,  with  a 
known  De/Ds  along  the  Sears  body  and  stress  gradients  taken  from  the 
Sears  body,  integration  was  allowed  to  proceed  from  the  stagnation 
poi nt. 

More  specifically,  the  same  local  similarity  criteria  were  employed 
to  define  a  point  on  a  Sears  body.  The  right-hand  side  of  Eq.  (31)  and 
the  pressure  gradients  were  then  calculated.  Moreover,  from  the  exact 
solution  of  the  Sears  problem  a  relationship  between  3xx/3z  and  3x^/3z 
was  found. 

!r  ■  u.k(2)  ♦  if  **«;<* »  <34> 


These  two  equations  were  then  solved  for  the  stress  gradients  assuming  a 
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constant  small  value  for  the  quantity  pq  in  the  denominator  of  Eq. 

(31).  The  stress  gradients  were  then  used  in  the  integration  of  Eq. 

(31)  over  the  prolate  spheroid,  assuming  the  same  small  value  for  the 
2 

quantity  pq  .  It  was  found  that  the  results  were  not  sensitive  to  the 

2 

particular  value  of  pq  .  Results  obtained  with  this  method  are  plotted 
in  Fig.  5.  It  is  now  interesting  to  observe  that  the  skin  friction 
lines  converge  to  an  envelope  on  both  sides  of  the  line  of  open  separa¬ 
tion. 

Yet  another  method  was  discovered  by  one  of  us  (C.C.)  who  made  the 

observation  that  on  all  problems  considered  by  a  variety  of  methods  the 

inviscid  velocity  is  constant  along  the  line  of  separation  to  within  an 

accuracy  of  a  few  percent.  This  value  can  be  estimated  from  a  Sears 

solution,  although  a  more  accurate  value  could  be  obtained  by  a  more 

accurate  interactive  method.  For  a  unit  oncoming  flow,  the  potential 
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velocity  on  a  Sears  body  takes  the  value  3.71 7sin  a  +  cos  a,  where  a  is 
the  yaw  angle.  Searching  for  the  locus  of  points  where  the  potential 
velocity  takes  this  value  on  a  prolate  spheroid  one  falls  to  within  3% 
of  the  separation  line  calculated  by  sophisticated  methods.  To  improve 
on  this  idea,  it  was  attempted  to  include  a  correction  term: 

Vgep  =  3.717sin2<x  +  cos2<x  +  f(p  a)  (36) 

where  f  may  depend  on  the  axes  ratio  and  the  yaw  angle. 

A  good  choice  for  the  function  f  can  be  made  by  considering  the 
limiting  values  of  a  and  a/b.  Clearly  for  a  ■  0°  and  a  ■  90°  the  prob¬ 
lem  becomes  truly  axisymmetric  and  truly  two-dimensional,  respectively. 


A  reasonable  value  for  a  =  0  is  the  velocity  at  separation  over  a 
sphere*®:  V  =  sin!09.6)  =  1  .9968  and  thus  by  comparison  to  Eq. 

(28),  f ( 1 .0 )  =  sin  (41.72)  =  0.99647.  A  reasonable  assumption  for 
f  then  could  be 


f(f,  a) 


tl\2  sin 
'2' 


(41. 72 -a 
9(a/b) 


The  function  g  should  have  the  following  limits 

g(a/b)  ♦  •  as  a/b  -►  ® 
g(a/b)  ♦  1  as  a/b  1 

since  a/b  »  1  and  a/b  =  »  correspond  to  a  sphere  and  an  infinite  cylin¬ 
der,  respectively.  A  comparison  with  detailed  earlier  calculations  for 
prolate  spheroids  indicates  that  for  values  of  a/b  equal  to  1,4  and  6, 
the  function  g(a/b)  should  take  the  values  0.89,  3.78  and  8.75,  respec¬ 
tively.  A  function  that  would  fit  through  such  data  is 

g(a/b)  «  0.1031 (a/b)2*43  +  0.7869  (38) 

Formula  (36)  with  f  and  g  given  by  Eqs.  (37)  and  (38)  was  used  to 
obtain  separation  lines  on  prolate  spheroids.  The  results  of  this 
method  and  the  method  of  local  similarity  described  earlier  were  com¬ 
pared  against  results  obtained  by  other  methods.  Results  for  axes 
ratios  of  a/b  *  2,4  and  6  and  various  angles  of  attack  are  plotted  in 
Figs.  6,  7  and  8,  respectively.  The  agreement  of  the  present  crude 

O  IQ 

method  with  the  methods  of  Stock  and  Geissler  is  surprisingly  good. 

Comparison  with  the  work  of  Refs.  3  and  6  is  given  in  Figs.  9,  10 
and  11.  In  Figs.  10  and  11  our  estimates  of  the  separation  line  as 
obtained  from  our  flow  visualizations  are  also  included.  For  a  =  6°  no 
experimental  data  is  available.  However,  for  a  *  5°,  as  will  be  dis¬ 
cussed  later,  there  is  no  experimental  evidence  of  closed  separation. 
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3.  Flow  Visualization 

Experiments  were  conducted  on  a  prolate  spheroid  model  with  axes 
ratio  b/a  =  1/4.  The  model  was  machined  out  of  aluminum  and  was  mounted 
by  a  strut  from  the  one  end.  Dye  ports  were  drilled  at  strategically 
chosen  points.  Colored  dyes  were  supplied  by  small  flexible  hoses  which 
were  fed  and  controlled  by  hypodermic  needles.  A  novel  method  was 
introduced  in  the  present  studies.  Thin  hypodermic  needles  were  mounted 
on  the  skin  of  the  body,  allowing  the  dyes  to  be  released  at  the  edge  of 
the  boundary  layer.  In  this  way,  skin  friction  lines  and  inviscid 
streamlines  emanating  from  the  same  station  in  space  could  be  compared. 
The  tests  were  conducted  in  a  water  tunnel,  and  still  color  photographs 
were  taken  from  views  corresponding  to  inclinations  <t>  *  0°,  90°,  120° 
and  180°.  These  correspond  to  windward,  side,  oblique  and  leeward  views, 
respectively. 

Experiments  were  conducted  for  angles  of  attack  a  *  3°,  6°,  10°, 
15°,  20°  and  30°.  The  results  confirmed  the  fact  that  for  very  low  and 
very  high  angles  of  attack,  here  for  a  3  3°  and  for  a  >  30°,  the  line  of 
separation  is  closed  around  the  body.  In  this  case,  two  critical  points 
can  be  identified  on  the  line  of  separation  and  the  general  pattern  fits 
nicely  one  of  the  models  described  by  Maskell  (see  discussion  in  Ref. 

3),  namely  the  bubble  model.  This  in  the  more  common  terminology  is  the 
case  of  closed  separation. 

In  Figs.  12-15  we  display  flow  visualizations  of  different  views  of 
the  flow  about  a  prolate  spheroid  at  angles  of  attack  of  3°,  6°,  10°  and 
20°,  respectively.  In  Fig,  12,  a  well-behaved  closed  separation  pattern 
is  identified  in  agreement  with  all  analytical  predictions. 


15 


twsi 


For  angles  of  attack  6°  <  a  <  30°  the  evidence  clearly  indicates 
that  there  are  two  lines  of  open  separation.  However,  details  of  the 
flow  uncovered  for  the  first  time  here  shed  light  to  controversial 
interpretations  of  earlier  numerical  results.  For  this  range  of  angles 
of  attack,  it  appears  that  the  basic  skin-friction  line  pattern  is  qual¬ 
itatively  the  same.  Three  distinct  regions  difficult  to  be  accessed  by 
the  free  stream  become  quickly  apparent  and  may  be  erroneously  inter¬ 
preted  as  regions  of  separated  flow.  The  first  region  is  covering  the 
lee  side  and  appears  dark  in  the  top  view  of  Fig.  14.  It  should  be 
emphasized  that  all  the  darkened  area  in  this  figure  is  due  to  dye  that 
has  been  emitted  at  only  one  port  in  the  front  part  of  the  model.  In 
the  same  figure  the  two  lighted  areas  on  the  two  sides  in  the  aft  of  the 
body  are  strongly  suggesting  the  existence  of  two  separated  regions. 

The  pattern  could  be  sketched  approximately  as  in  Fig.  16.  The  two  side 
lobes  would  then  appear  to  have  the  shape  described  by  Cebeci  et  al6. 
What  has  been  known  for  a  long  time  is  that,  in  fact,  the  flow  remains 
attached  on  the  leeward  plane  of  symmetry  and  in  its  neighborhood.  This 
is  corroborated  by  the  fact  that  the  dye  stays  on  the  wall  until  it 
reaches  the  small  mount  at  the  aft  of  the  body.  However,  our  evidence 
indicates  that  what  appear  as  side-separation  lobes  are  actually  regions 
of  incredibly  large  cross  flows.  This  is  shown  more  clearly  in  Fig.  17. 
In  this  figure  the  inviscid  streamlines  as  well  are  marked  by  dyes 
released  at  the  edge  of  the  boundary  layer. 

A  very  large  number  of  flow  visualizations  obtained  point  very 
clearly  to  a  pattern  of  viscous  streamlines  shown  in  Figs.  18  and  19 
for  o  *  10°  and  20°,  respectively.  These  patterns  have  been  drawn  care¬ 
fully  by  an  artist  who  traced  the  edges  of  the  dye  lines  from  a  large 
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number  of  flow  visualizations.  We  are  presently  in  the  process  of  using 
a  pattern-recognition  computer  technique  to  store  all  this  information 
in  digital  form  and  eventually  plot  the  equivalent  of  Figs.  18  and  19  by 
purely  mechanical  means. 

Our  evidence  as  shown  in  Figs.  13,  14,  15  and  the  interpretations 
shown  in  Figs.  18  and  19  indicates  that  for  these  angles  of  attack,  the 
separation  line  is  of  the  open  type.  In  fact,  the  dyes  lift  off  from 
the  body  only  along  the  envelope  of  the  skin  friction  line.  This  has 
been  observed  by  the  present  investigators,  but  so  far  it  appeared 
difficult  to  photograph  except  in  Fig.  17.  However,  the  inviscid 
streamlines  of  Figs.  15  and  17  also  indicate  that  the  inviscid 
streamline  in  the  neighborhood  of  the  separation  line  lift  off  by  a 
rather  large  angle  with  respect  to  the  body. 

In  view  of  this  evidence,  we  feel  that  Cebeci  et  al^  simply  reached 
or  rather  came  very  close  to  the  point  of  zero  axial  shear  stress  when 
they  were  forced  to  abandon  their  integration  on  the  leeside  of  the 
body.  This  possibility  was  already  identified  in  their  presentation. 
Apparently,  the  locus  of  such  points  may  be  impossible  to  pass  over  with 
conventional  numerical  schemes  or,  in  the  terminology  of  Ref.  6,  such 
flows  are  incalculable.  Evidence  from  the  flow  visualizations  presented 
here  as  well  as  the  very  large  number  of  data  which  are  not  included, 
indicate  that  in  this  region  the  crossflow  angles  may  reach  values  of 
the  order  of  150°.  It  is  our  strong  conviction  that  in  the  case  studied 
in  Ref.  6,  the  skin-friction  lines  have  the  qualitative  pattern  shown  in 
Fig.  20. 


4.  Velocity  Measurements 

In  this  section  we  describe  a  method  of  measuring  velocity  profiles 
in  three-dimensional  boundary  layers.  It  is  emphasized  here  once  more, 
that  to  the  knowledge  of  the  present  authors,  Laser-Doppler  Velocimetry 
has  not  yet  been  used  in  measuring  three-dimensional  flows.  Moreover, 
with  the  exception  of  Meier's  group,  12-15  wh0  velocities  with 

3-hole  probes,  no  one  has  measured  velocity  profiles  on  a  prolate 
spheroid.  The  present  results  should  therefore  provide  useful  informa¬ 
tion  to  analyticians  who  encounter  great  difficulties  in  the  lee  side  of 
a  prolate  at  an  angle  of  attack. 

A  special  Laser-Ooppler  Velocimeter  was  designed  and  constructed 
for  the  present  project.  The  basic  consideration  for  such  a  design  is 
that  at  least  two  components  of  velocity  are  needed  in  three-dimen¬ 
sional-boundary  layers  to  generate  information  in  the  streamwise  and  the 
crossflow  direction.  Moreover,  these  components  must  be  parallel  to  the 
plane  tangent  to  the  body  surface  at  the  point  of  interest.  It  is  also 
necessary  that  the  measuring  volume  be  allowed  to  approach  the  wall 
without  blockage  of  any  of  the  beams  by  the  model.  To  meet  these 
requirements,  three  laser  beams  were  generated  in  the  same  plane,  as 
shown  in  Fig.  21.  With  proper  polarization  of  these  beams  it  is  pos¬ 
sible  to  receive  scattering  signals  from  interference  of  beams  1  and  2, 
and  2  and  3.  In  this  way  the  components  of  the  velocity  in  the  direc¬ 
tions  of  the  bisectors  can  be  obtained.  A  computer  program  then  allows 
the  determination  of  Cartesian  components. 


A  special  optical  bench  was  constructed.  The  bench  was  mounted  on 
traversing  towers  which  permit  rotation  about  the  optical  axis  of  the 
LDV  as  well  as  translation  in  the  horizontal  and  vertical  direction  as 


shown  in  Fig.  22.  The  train  of  sending  optics  consists  of  two  beam 
splitters.  The  first  is  combined  with  a  polarization  rotator  to  allow 
control  of  the  power  of  the  two  beams.  Bragg  cells  may  be  used  and,  in 
fact,  will  be  necessary  to  allow  measurements  in  regions  of  reversing 
flows.  The  beams  are  then  all  brought  in  one  plane.  The  outer  two 
beams  are  displaced  outward  and  then  turned  back  to  meet  at  the  measur¬ 
ing  volume. 

A  key  feature  of  the  present  system  is  that  traversing  is  permitted 
perpendicular  to  the  local  tangent  plane  to  the  body.  This  is  the 
natural  boundary  layer  coordinate,  normal  to  the  surface  of  the  body. 
Moreover,  the  two  components  parallel  to  the  wall  surface  can  be  mea¬ 
sured  all  the  way  down  to  the  wall.  This  is  possible  because  all  three 
beams  are  in  one  plane.  A  stepping  motor  has  been  connected  to  allow 
accurate  traversing  across  the  boundary  layer. 

A  crucial  step  in  Laser-Doppler  measurements  is  the  determination 
of  the  location  of  the  measuring  volume.  In  the  present  experiments 
this  was  accomplished  as  follows.  The  location  of  the  model  in  the  test 
section  was  first  determined  optically  by  virtue  of  a  telescope.  The 
space  coordinates  of  the  two  points  of  the  spheroid  surface  along  its 
axis  of  symmetry  were  recorded.  The  optical  bench  was  then  displaced 
and  rotated  until  the  plane  of  the  beams  became  tangent  to  the  model  and 
the  point  of  the  beam  intersection  was  brought  as  close  as  possible  to 
the  model  surface.  The  locations  of  the  intersections  of  the  three 
beams  with  the  two  side  walls  of  the  tunnel  were  then  determined  and 
stored.  A  computer  program  was  prepared  to  estimate  the  location  of  the 
measuring  volume  on  the  model  by  reconstructing  the  model  and  the  beams. 
An  example  of  this  geometry  is  shown  in  Fig.  23.  The  method  allows  the 


operator  to  position  the  system  as  close  as  possible  to  the  desired 
point.  The  final  position  of  the  measuring  volume  on  the  model  is  then 
recorded.  Traversing  across  the  boundary  layer  is  performed  perpendic¬ 
ular  to  the  model  body  in  steps  controlled  by  the  stepping  motor. 

It  was  estimated  that  the  error  of  positioning  the  measuring  volume 
on  the  body  in  ci rcumferentai 1  or  axial  location  should  be  less  than  1% 
of  the  length  of  the  model.  Errors  in  vertical  displacements  should  be 
less  than  0.1%  of  the  average  thickness  of  the  boundary  layer.  It 
should  be  emphasized  that  the  actual  location  of  the  wall,  namely  the 
origin  of  the  z-axis  was  determined  by  extrapolation  after  all  the  data 
were  obtained. 

A  third  channel  consisting  of  a  second  laser,  beam  splitter,  col¬ 
lecting  and  receiving  optics  has  been  added  to  the  system  for  measure¬ 
ment  of  the  third  component  of  the  velocity.  Four  beams  thus  converge 
at  the  measuring  volume  as  shown  in  Fig.  24.  These  are  not  orthogonal 
to  each  other,  but  Cartesian  components  can  be  calculated  on  line  by  the 
laboratory  computer,  in  terms  rf  the  angles  shown  in  the  figure. 

An  automatic  system  has  been  designed  and  tested  for  collecting  and 
storing  all  three  components  for  as  many  points  in  the  boundary  layer  as 
necessary.  This  method  has  been  tested,  but  only  two  components  were 
obtained  and  reported  here.  With  speeds  as  low  as  3  to  10  cm/sec,  it 
was  possible  to  generate  boundary  layers  with  thicknesses  of  a  few 
millimeters  on  a  model  17  cm  long.  We  were  able  to  obtain  40  to  50  data 
points  in  such  boundary  layers,  but  600  or  700  points  can  be  easily 
obtained  by  cutting  the  traversing  step. 

The  experimental  setup  has  been  arranged  and  tested  together  with 
the  data  acquisition  system  on  a  small  water  tunnel  which  did  not  have 
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good  flow  quality.  It  was  planned  to  move  the  entire  experimental  rig 
and  reassemble  it  around  the  ESM  water  tunnel  for  obtaining  the  experi¬ 
mental  data  as  described  in  the  proposal.  Unfortunately,  a  sequence  of 
serious  problems  has  set  back  the  work  in  our  laboratory.  It  was  then 
decided  to  invest  time  and  effort  in  improving  the  quality  of  the  flow 
in  the  small  tunnel  in  order  to  complete  the  work  on  this  project  within 
the  alloted  time.  A  large  number  of  configurations  involving  combina¬ 
tions  of  honeycombs  and  screens  has  been  tested.  In  Fig.  25  we  show  the 
frequency  spectrum  of  the  velocity  RMS  signal.  We  were  very  successful 
in  reducing  drastically  the  turbulence  level  across  the  entire  scale, 
except  for  the  smallest  frequencies.  Apparently,  a  very  low  frequency 
is  unavoidable  and  is  due  to  the  finite  number  of  blades  of  the 
impeller.  We  then  replaced  the  driving  system  by  a  centrifugal  pump  and 
installed  a  number  of  foam  barriers  to  reduce  the  head  and  eliminate  the 
large  scale  vortices.  With  this  design  we  were  able  to  lower  the  over¬ 
all  turbulence  level  to  approximately  0.35%. 

The  heart  of  the  data  acquisition  system  was  a  very  inexpensive 
time  switch.  This  was  dictated  by  the  large  power  requirements  of  sole¬ 
noids  that  control  blinders  of  the  laser  beams.  Measurements  were 
obtained  sequentially  as  follows.  A  blinder  blocked  beam  1  and  there¬ 
fore  the  component  V23  (see  Fig.  21)  could  be  measured.  The  time  switch 
sent  a  message  to  the  computer  to  obtain  a  measurement.  Beam  3  was  then 
blocked,  and  the  computer  was  signaled  to  measure  and  store  the  value  of 
the  component  V^.  A  third  signal  from  the  time  switch  to  the  computer 
led  the  program  to  send  pulses  to  the  stepping  motor.  The  number  of 
pulses  is  preprogrammed  to  allow  fine  steps  near  the  wall. 


Laser  beam  2,  the  optical  axis  of  symmetry  of  the  Laser  Veloci meter 
was  always  perpendicular  to  the  direction  of  the  free  steam.  However, 
this  axis  was  perpendicular  to  the  ellipsoid  axis  for  only  one  position 
of  the  ellipsoid.  This  is  because  the  model  was  inclined  with  respect 
to  the  free  stream  by  the  angle  of  attack,  a,  but  was  rotated  about  an 
axis  parallel  to  the  free  stream  to  expose  all  surfaces  from  the  wind¬ 
ward  to  the  leeward  sid4  to  the  LDV  beams.  As  a  result,  to  obtain  velo¬ 
city  components,  it  was  necessary  to  measure  the  angle  4  between  the  x- 
axis  and  the  axis  of  the  LDV  system  as  shown  in  Fig.  21.  The  velocity 

components  u  and  v  (see  Fig.  1  for  definition)  are  then  given  in  terms 

of  the  measured  quantities  and  the  angles  defined  in  Fig.  21  by  virtue 
of  the  following  straightforward  relationships: 

V12  =  u  cosei2  +  v  sin0i2  (38) 

V23  =  u  cos913  +  v  sin  e13  (39) 

Errors  in  the  magnitude  of  the  velocity  on  the  order  of  0.5%  of  the 
mean  may  be  anticiapted.  However,  a  simple  analysis  indicates  that  for 
an  included  angle  of  0  =6°  (see  Fig.  21),  the  error  anticipated  in  the 
circumferential  direction  may  be  up  to  5%  of  the  mean  axial  velocity. 
Considering  the  fact  that  the  velocity  component  in  this  direction  may 
be  rather  low,  this  error  is  sometimes  considerable.  Clearly  a  larger 
angle  4  should  be  used  in  future  experiments. 

A  standard  cubic-spline  smoothing  routine  was  employed  to  generate 
velocity  profiles  from  the  discrete  experiment  data.  Comparison  of  the 
smoothed  profiles  and  actual  experimental  data  are  shown  in  Fig.  26.  In 
this  figure  it  is  also  apparent  that  some  data  are  missing  in  the 
immediate  vicinity  of  the  wall.  The  inability  to  obtain  data  in  this 
region  was  attributed  to  the  fact  that  the  size  of  the  measuring  volume 


was  not  very  small  compared  to  the  thickness  of  the  boundary  layer.  The 
situation  is  more  critical  for  axial  stations  near  the  nose  where  the 
thickness  of  the  boundary  layer  is  very  small.  In  such  regions  up  to 
10%  of  the  velocity  profile  are  void  of  experimental  data.  The  problem 
can  be  remedied  by  using  smaller-size  pinholes  in  the  photomultipliers. 

Data  were  obtained  in  this  way  at  five  axial  stations  and  five  cir- 
cumferentail  directions.  To  facilitate  comparison  with  the  numerical 
results  of  Ref.  20,  the  present  experimental  results  are  presented  in  an 
outlay  similar  to  the  one  included  in  this  reference.  The  present 
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experimental  results  and  the  numerical  results  of  Patel  &  Baek  are 
shown  in  Figs.  27  and  28,  respectively.  It  should  be  noticed  that  the 
position  of  the  experimental  profiles  does  not  coincide  with  integral 
values  of  coordinates.  This  is  due  to  difficulties  encountered  in  posi¬ 
tioning  the  measuring  volume  on  the  surface  of  the  model.  The  origin  of 
the  velocity  profile  in  this  figure  is  positioned  at  the  point  in  the 
x-<t»  plane  with  the  coordinates  of  the  z=0  position  of  the  measuring 
volume. 

on 

Patel  and  Baektu  found  very  abrupt  changes  in  the  thickness  of  the 
boundary  layer  and  the  shape  of  the  velocity  profile  in  the  neighborhood 
of  the  line  of  separation,  namely,  between  the  angles  0  *  130°  and  135°. 
The  velocity  profile  in  the  region  beyond  separation  shows  clear  indica¬ 
tions  of  an  inflectional  shape.  This  agrees  at  least  qualitatively  with 
our  experimental  data.  Indeed,  the  first  inflectional  profile  appears 
at  <j>  =  130°  but  not  until  x/a  =  0.4.  However,  the  abrupt  change  of  the 
thickness  of  the  boundary  layer  is  not  present.  Apparently,  the  bound¬ 
ary  layer  equations  are  an  insufficient  model  for  describing  exactly  the 


physical  situation  in  the  neighborhood  of  separation  much  like  the  two- 
dimensional  case. 

20 

The  careful  numerical  calculations  of  Patel  and  Baek  indicate  in 
the  forward  half  of  the  body  thicker  boundary  layers  closer  to  the 
windward  side  than  the  experiments.  Indeed,  the  same  trend  is  indicated 
by  the  calculations  of  Cebeci  et  al  ,6  as  shown  in  Figs.  9  and  10. 
However,  this  pattern  is  reversed  for  the  aft  portion  of  the  body. 

Both  components  of  the  velocity  in  the  axial  and  circumferential 
direction  are  plotted  in  Figs.  29-32  for  four  axial  stations  along  the 
model.  The  reader  is  cautioned  to  the  fact  that  in  these  figures  five 
different  scales  of  measurement  are  being  employed.  The  position  x,<j>  of 
the  point  on  the  ellipsoid  where  the  data  were  obtained  is  the  origin  of 
the  profile  axis.  The  detailed  markings  on  the  orthogonal  system  shown 
in  each  figure  provide  the  velocity  scales  and  the  scale  in  the  vertical 
axis.  The  two  velocity  scales  parallel  to  the  x  and  $  directions, 
respectively,  should  not  be  confused  with  the  corresponding  length 
scales  in  these  directions.  Both  velocity  scales  are  marked  in  tenths 
of  the  free  stream  velocity,  whereas  the  unit  length  in  the  z  axis  cor¬ 
responds  to  1/10  in.  or  0.85  mm. 

Finally,  a  collection  of  all  the  experimental  data  obtained  is  pre¬ 
sented  in  a  three-dimensional  outlay  in  Fig.  33.  This  figure  is  a  com¬ 
posite  of  Figs.  29-32  and,  therefore,  the  scales  are  identical  to  those 
employed  previously. 


5.  Conclusions 

The  analytical  methods  presented  here  have  been  proven  to  provide 
remarkable  accuracy  with  a  minimum  of  computer  time  and  space.  The 
experimental  results,  on  the  other  hand,  seem  to  reconcile  the  points  of 
view  of  analyticians  about  closed  and  open  separation.  It  appears  that 
for  a  good  range  of  angles  of  attack  there  exists  open  separation  which 
is,  in  fact,  compatible  with  boundary-layer  studies.  The  present  work 
generated  a  large  number  of  flow  visualizations  which  conclusively  prove 
this  point. 

The  LDV  measurements  were  obtained  at  the  five  stations  in  the  for¬ 
ward  portion  of  the  ellipsoid.  However,  it  appears  necessary  to  obtain 
more  data  further  downstream  and,  in  fact,  more  detail  is  necessary  in 
the  neighborhood  of  the  origin  of  open  separation. 

The  work  described  here  is  the  result  of  one  year's  work  according 
to  the  plan  described  in  our  contract  with  the  Navy.  However,  the  find¬ 
ings  clearly  indicate  the  need  for  more  work  to  combine  and  interpret 
correctly  flow  visualizations  and  LDV  measurements  in  the  areas  of  great 
interest. 

6.  Acknowlegment 

The  support  of  the  Naval  Air  Systems  Command  (310-D)  under  contract 
No.  N00167-82-K-0085  is  greatly  appreciated.  The  authors  are  also 
indebted  to  the  monitor  of  the  work.  Dr.  T.  C.  Tai  (DTNSRDC-1606) ,  for 
his  personal  and  continuous  encouragement  and  his  valuable  technical 


assistance 


References 


1.  Blottner,  F.  G.,  and  Ellis,  M.  A.,  "Finite-difference  solution  of 
the  incompressible  three-dimensional  boundary-layer  equations  for  a 
blunt  body".  Computers  &  Fluids,  Vol.  1,  pp.  133-158,  1973. 

2.  Geissler,  W.,  "Three-dimensional  laminar  boundary  layer  over  a  body 
of  revolution  at  incidence  and  with  separation",  AIAA  J.,  Tech. 
Note,  Vol.  12,  No.  12,  pp.  1743-1745,  1974. 

3.  Wang,  K.  C.,  "Boundary  layer  over  a  blunt  body  at  low  incidence 

with  ci rcumferencial  reversed  flow",  J.  Fluid  Mech.,  Vol.  72,  part 
1,  pp.  49-65,  1975.  “  ~  ' 

4.  Hall,  M.  G.,  "A  numerical  method  for  calculating  steady  three- 
dimensional  laminar  boundary  layers",  RAE  TR  G7145,  1967. 

5.  Patel,  V.  C.,  and  Choi,  D.  H.,  "Calculation  of  three-dimensional 
laminar  and  turbulent  boundary  layers  on  bodies  of  revolution  at 
incidence",  2nd  Symposium  on  Turbulent  Shear  Flows,  pp.  15.14- 
15.24,  Imperial  College,  London,  1979. 

6.  Cebeci ,  T.,  Khattab,  A.  K.,  and  Stewartson,  K.,  "Three-dimensional 
laminar  boundary  layers  and  the  OK  of  accessibility",  J.  Fluid 
Mech.,  Vol.  107,  pp.  57-87,  1981. 

7.  Ragab,  S.  A.,  "A  method  for  the  calculation  of  three-dimensional 
boundary  layers  with  circumferential  reversed  flow  on  bodies", 
AIAA-82-1023,  1982. 

8.  Stock,  H.  W.,  "Computation  of  the  boundary  layer  and  separation 
lines  on  inclined  ellipsoids  and  of  separated  flows  on  infinite 
swept  wings",  AIAA-80-1442,  1980. 

9.  Tai,  T.  C.,  "Determination  of  three-dimensional  flow  separation  by 
a  streamline  method",  AIAA  J.,  Vol.  19,  No.  10,  Oct.  1981,  pp  1264- 
1271. 

10.  Tai,  T.  C.,  "Effect  of  cross  flow  on  the  vortex-layer-type  3-dimen¬ 
sional  flow  separation",  in  Three-Dimensional  Turbulent  Boundary 
Layers,  eds.  H.  H.  Fernholz  and  E.  Krause,  pp.  117-125,  1^82. 

11.  Han,  T.  and  Patel,  V.  C.,  "Flow  separation  on  a  spheroid  at  inci¬ 
dence",  J.  Fluid  Mech.,  Vol.  92,  Part  4,  pp.  643-657,  1979. 

12.  Kreplin,  H.-P.,  Volmers,  H.,  and  Meier,  H.  V.,  "Experimental  deter¬ 
mination  of  wall  shear  stress  vector  on  an  inclined  prolate 
spheroid",  AFFDL-TR-80-3088,  1980. 

13.  Meier,  H.  U.,  and  Kreplin,  H.-P.,  "Experimental  investigation  of 
the  transition  and  separation  phenomena  on  a  body  of  revolution", 

Z.  Flugwiss.  Weltraumforsch.  4,  Heft2,  pp.  65-71,  1980. 


14.  Meier,  H.  U.,  adn  Kreplin,  H.-P.,  "Experimental  study  of  boundary 
layer  velocity  profiles  on  a  prolate  spheroid  at  low  incidence  in 
the  cross  section  x/L  =  0.64",  AFFDL-TR-80-3088,  pp.  169-189,  1980 

15.  Meier,  H.  U.,  Kreplin,  H.-P.,  and  Vollmers,  H.,  "Velocity  distribu 
tions  in  3-0  boundary  layers  and  vortex  flows  developing  on  an 
inclined  prolate  spheroid",  DFVLR-AVA-REPORT  IB  222  81  CP  1,  1981. 

16.  Sears,  W.  R.,  "The  boundary  layers  of  yawed  cylinders",  J.  Aero. 
Sci.,  Vol .  15,  pp.  49-52,  1942. 

17.  Band,  E.  6.  U.  and  Payne,  P.  R.,  "The  Pressure  Distribution  on  the 
Surface  of  an  Ellipsoid  in  Inviscid  Flow",  Aeronautical  Quarterly. 
pp.  70-84,  February  1979. 

18.  N.  Frossling,  Verdunstung,  Warmeubergang  und  Geschwindi gkeitsver- 
teilung  bei  zwei -dimensionaler  und  rotationssymmetrischer 
laminarer,  Greuzschichtstromung.  Lunds.  Univ.  Arsskr.  N.F.  Avd.  2, 
35,  No.  4,  1940. 

19.  Geiss.ler,  W.,  "Berechnung  der  laminaren  Ablosung  an  unsymmetrisch 
umstromten  Rotationskorpern  mittels  Differenzenverfahren,  DFVLR-1B 
pp.  251-273,  A-31 ,  1973. 

20.  Patel,  V.  C.  and  Baek,  J.  H.,  "Calculation  of  boundary  layers  and 
separation  on  a  spheroid  at  incidence",  in  Numerical  &  Physical 
Aspects  of  Aerodynamic  Flows,  ed.,  T.  Cebeci ,  Cal.  State  Univer¬ 
sity,  1983. 

21.  Ramaprian,  B.  R.,  Patel,  V.  C.  and  Choi,  D.  H.,  "Mean-flow  measure 
ments  in  the  three-dimensional  boundary  layer  over  a  body  of  revo¬ 
lution  at  incidence",  J.  Fluid  Mech.,  Vol.  103,  pp.  479-504,  1981. 


Lines  of  separation  for  a  prolate  spheroid  with  axes  ratios  a/b 

s  6  and  angles  of  attack  a  =  5°,  10°  and  15°.  _ ,  Stock8, 

. »  Geissleriy;  □  ,  present  method. 


7  Lines  of  separation  for  a  prolate  spheroid  with  axes  ratios  a/b 

=  4  and  aggies  of  attack  a  =  15°  and  30°.  - ,  Stockd;  - 

Geissler19; n,  present  method. 


Cf  s0— -\  Accessible  line* 

-  e  d\ 


External  Streamlines 


r-Separation  line 
( Ref.  6 ) 

■Separation  line 
present 

_  measurements 


(0  ' 


Fig.  10  Lines  of  separation  and  external  streamlines  for  a  prolate 

spheroid  of  a/b  =  6  at  an  angle  of  attack  a  =  15°  as  reported 
in  Ref.  6;  Q,  present  approximate  method;  -o-o-o-o;  present 
experimental  data. 
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Fig.  11  Lines  of  separation  and  external  streamlines  for  a  prolate 

spheroid  of  a/b  =  4  at  an  angle  of  attack  a  =  30°  as  reported 
in  Ref.  6;  Q  ,  present  approximate  method;  o-o-o-o-o;  present 
experimental  data. 


Fig.  12  Skin  friction  lines  for  a  -  3°.  (a)  Leeward  view,  (b)  View 

from  $  =  120°,  (c)  Side  view. 


Fig. 21  Schematic  representation  of  the  crossing  of  3  beams.  All  3  beams 
are  in  the  same  plane.  Interference  between  beam  1  and  2  generates 
the  velocity  component  V,?  and  interference  between  2  and  3  generates 
the  component  V„-,.  * 
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23  Example  of  laser  beam  locations  with  respect  to  the  body 
and  output  of  program  for  the  calculation  of  measuring 


position. 
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Fig.  25  RMS  spectrum  of  the  velocity  fluctuation.  If  the  low  frequen¬ 
cies  (less  than  1  Hz)  are  eliminated,  then  the  turbulence  level 
Is  less  than  0.1%.  The  spikes  In  the  spectrum  correspond  to 
the  60  Hz  and  its  multiples  and  are  due  to  the  power  supply. 
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The  u-component  velocity  profile.  Horizontal 
scales  are  marked  in  tenths  of  free  steam  vel 
ments  of  1/10"  (and  0.85  mm)  respectively. 


The  u-  and  v-velocity  components  at  the  first  station  of  mea 
surement. 
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